Exact results for classical Casimir interactions: 
Dirichlet and Drude model in the sphere-sphere and sphere-plane geometry 
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Analytic expressions that describe Casimir interactions over the entire range of separations have 
been limited to planar surfaces. Here we derive analytic expressions for the classical or high- 
temperature limit of Casimir interactions between two spheres (interior and exterior configurations), 
including the sphere-plane geometry as a special case, using bispherical coordinates. We consider 
both Dirichlet boundary conditions and metallic boundary conditions described by the Drude model. 
At short distances, closed-form expansions are derived from the exact result, displaying an intricate 
structure of deviations from the commonly employed proximity force approximation. 



PACS numbers: 12.20.-m, 03.70.+k, 42. 25. Fx 

The collective action of fluctuation "van der Waals" 
forces between individual atoms leads to forces between 
macroscopic surfaces p] [2] • The fluctuations can be of 
thermal [5] and/or quantum [3] nature. An important 
characteristic of these forces is their non-additivity which 
complicates the study of Casimir interactions in the ther- 
modynamic limit involving a macroscopic number of par- 
ticles. An important exception is the interaction of pla- 
nar, parallel surfaces where symmetry allows for an ex- 
act solution, the so-called Lifshitz formula [5 , which is a 
landmark in the physics of fluctuation forces. More com- 
plicated shapes have been studied at sufficiently short 
surface separations by the proximity force approxima- 
tion (PFA) that is based on the Lifshitz formula, treat- 
ing curvature by summing over planar surface elements 
[5J. The thermal Casimir interaction between spherical 
particles in a critical fluid has been computed using con- 
formal invariance for nearly touching and widely sepa- 
rated particles [7] . More recently, a number of analytical 
and numerical approaches have been developed to deal 
with the non-additivity, most notably implementations 
of concepts from scattering theory 8 TO . However, it 
remains unclear to what extend and precision these ap- 
proaches can handle the practically important limit of 
short distances. Attempts to obtain analytic predictions 
beyond the PFA have been limited to compute first order 
corrections from a gradient expansion |11H13) . 

The difficulty to compute deviations from the simple 
PFA can be traced back to the long-range nature of fluc- 
tuation forces which is poorly treated by the PFA. With 
increasing spatial dimensions, fluctuations decay more 
strongly, and the PFA is expected to become more re- 
liable. The situation resembles the situation in statis- 
tical mechanics when phase transitions are studied by 
mean field theory, treating fluctuations poorly. In this 
context, the exact solution of the two-dimensional Ising 



model [13] helped understanding phase transitions and 
inspired several developments in the theory of critical 
phenomena and related fields. This demonstrates the 
importance of exact solutions for a better understand- 
ing of complex phenomena. In particular, several forms 
for non-additivity induced deviations from the PFA have 
been hypothesized before, depending on boundary con- 
ditions and temperature [HI 1151 H6j . Here we derive an 
exact solution for the interaction between spheres of dif- 
ferent radii that displays a much richer structure for the 
Casimir interaction in the high temperature limit than 
previously assumed. 

In this Letter, we obtain analytic expressions for the 
high-temperature limit of the Casimir interactions be- 
tween two spheres (interior and exterior configurations), 
including the sphere-plane geometry as a special case. 
We consider both a scalar field with Dirichlet boundary 
conditions and metallic boundary conditions described 
by the Drude model. The key tool of our approach 
is a bispherical coordinate system [17j . which allows 
for a simple solution of the scattering problem of two 
spheres in the static limit (Laplace equation) . The high- 
temperature sphere-plate problem has been analyzed be- 
fore in the large distane limit [18] and by large-scale nu- 
merics at short distances, including up to 5000 partial 
wave orders [T5]. Our exact solution is universal, i.e., 
material independent and displays full agreement with 
these two limiting cases. We derive from the exact result 
a closed-form short-distance expansions that reveals an 
intricate structure of deviations from the PFA. It turns 
out that for Dirichlet boundary conditions, the classical 
Casimir force can be expanded as a Laurent series for 
small surface-to-surface separations L. For Drude metal- 
lic boundary conditions, the structure of the corrections 
to PFA is more complicated than suspected before [15], 
as it involves a double logarithm of L, as well as pow- 
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ers of InL multiplied by powers of L. In both cases the 
leading correction to the PFA energy is represented by a 
term proportional to InL, with a common coefficient that 
had been computed earlier by using a recently pro- 
posed gradient expansion of the Casimir energy [T2J [T3] . 
However, for experimentally accessible separations, the 
interaction of Drude metals is dominated entirely by the 
double logarithmic term. 




FIG. 1: Geometry of two spheres and sphere-plate. Shown are 
the centers (C, C') of the spheres and their foci (F, F'). The 
dashed curves correspond to curves of constant bispherical 
coordinate /x with fi = =p7r, =F7r/2, =F7r/10, =F7r/20 with the 
upper (lower) sign for z < 0(z > 0), beginning with the 
smallest sphere. 
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Geometry, coordinates and eigenf unctions 
We consider two spheres of radii Ri and 
with surface-to-surface distance L. This 
ometry is conveniently parametrized in 
spherical coordinates (/x, 77, </?) pT], defined 
(x, y, z) — a(sin r\ cos if, sin r\ sin ip, sinh li) / (cosh /x— cos rj) 
where z = ±a are the foci of the two spheres defined by 
(X = fi% > and fx = fi 2 < 0, respectively, see Fig. [I] The 
spheres have radii R\ = a/sinh/xi and R 2 = — a/ sinh /i 2 
and distances L\ — a coth [i\ and L2 = —a coth \x 2 from 
the origin. The center-to-center distance is d = L\ + Li- 
lt is useful to express the fi a in terms of the natural ge- 
ometrical scales, fix — arccoshAi, \x 2 — — arccoshA2 with 
X a = [L 2 + 2(L + R a ){Ri + R 2 )]/[2R a (L + R x + R 2 )]. 
We are interested in the classical Casimir inter- 
action which is determined by the scattering of 
the spheres at zero frequency which is determined 
by the Laplace equation. The latter is separable 
in bispherical coordinates and its Green's func- 



tion can be expanded as Go(/i,?7, <p; (J.' ,?/ , f') = 

E to EL-i f> ^)C* , wh ere /i <(>) 

is the smaller (larger) of /1 and // and the regular and 
outgoing eigenfunctions are 



, rog/ out 



I cosh 



fi — cos 77 



a(2l + 1) 



Y lm (v,v)e ±il+1 ^ (1) 



for I > 0, m = — I, . . . , I. Scattering solutions for various 
boundary conditions can be expanded in these eigenfunc- 
tions, leading to the zero-frequency (static) scattering 
amplitude t/"|, / of sphere a by which we denote the 
matrix elements of the operator in the bispherical 
basis. The Casimir energy is expressed in the scattering 
approach [5] in terms of the scattering amplitudes and 
translation operators that translate the scattering solu- 
tion from the coordinate of one object to the one of the 
other object. Here, for both spheres the scattering am- 
plitude is expressed in the same bispherical coordinate 
system and hence the translation operators become the 
identity operators, yielding for the classical Casimir en- 
ergy (or th order Matsubara term) , 



E 



lndet[l-iV] with N = f (1) f (2) . (2) 



It is important to notice that the energy depends only 
on the equivalence class [[TV]] formed by all matrices 
that represent the operator N, i.e., two elements iV and 
MNM~ X of the class differ by an (invertible) similarity 
transformation M . In the following we denote by [[. . .]] 
the equivalence class of the matrix enclosed by the brack- 
ets. 

Dirichlet boundary conditions — The scattering ampli- 
tudes for Dirichlet boundary conditions tplp = fi a ) — 
on the spheres follows from Eq. ([IJ 8 . They assume 
the simple diagonal and m-independent form Ts~*\, , — 
— exp(=F(2Z + l)fx a )Su'S mm ' with — for a = 1 and + for 
a = 2. Hence, the determinant in Eq. ^ can be easily 
computed and we get for the classical Casimir energy of 
two Dirichlet spheres the exact result 



5])(2Z + l)ln 



2i + 11 



(3) 



which depends only via the single variable Z = [X% + 
a/ X\ — 1] _1 [A 2 + \/X\ — on the geometrical scales. 
For \Z\ < 1, including the range of physically mean- 
ingful parameters, the energy E^ D ) is an analytic func- 
tion [5T]. Large distances L correspond to small Z with 
Z = R1R2/L 2 + (D(L~ 3 ) and small distances to Z close 
to unity with Z = l- ^2L/R X + 2L/R 2 + 0{L). 

We note that the configuration of two spheres with 
one (i?i) inside the other (i? 2 ) [IH US can be com- 
puted exactly by the same method. The classical en- 
ergy is given again by Eq. Q but now with Z = 
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[Aa + v / A| r l]/[Ai + nAF 7 !] with Ai = [-L 2 + 2{L + 
Ri)(-Ri + R 2 )]/[2R 1 {-L - Rx + fl 2 )] > A 2 = [i 2 + 
2(-L + i? 2 )(-i? 1 +i? 2 )][2i? 2 (-L-i? 1 + i? 2 )] where L is 
the closest surface-to-surface separation of the spheres so 
that their centers have a distance d — i? 2 — Ri — L > 0. 
For two concentric spheres one has Z = i?i/i? 2 . 

Drude model — In general, a Drude metal has to be 
described within electromagnetic scattering theory. How- 
ever, in the classical limit, it can be shown that only the 
electric or TM modes contribute to the Casimir energy 
|18j . We shall use this observation below to map the in- 
teraction of Drude metals to a scalar field problem that is 
similar to the one with Dirichlet conditions. To simplify 
notations, we focus here on the sphere-plate geometry 
(with the sphere described by ^ = /ii > 0). When the 
operator N that yields the classical energy for the Drude 
model is expressed in a spherical wave basis with C as 
the origin (see Fig. [T]) , it has the form [18] 



N = 



(i + iy. 



(I + m)\{V - m)\ 



{z + z- 1 ) 



-l-l'-i 



(4) 



where we used that R/(L + R) = 2/(Z +1/Z) for a plate 
and a sphere of radius R. The same matrix is obtained for 
a scalar field with Dirichlet conditions but with one im- 
portant difference: While in the Dirichlet case monopoles 
with I = are included, in the Drude case monopoles are 
excluded since there are no corresponding spherical vec- 
tor waves. Physically, this feature is a consequence of 
the fact that an isolated compact object has zero total 
charge, i.e., no monopole. The next important observa- 
tion is that a given spherical multipole of order (l,m) 
decomposes into infinitely many bispherical multipoles 
(?', m) with m conserved and hence in the bispherical ba- 
sis the difference between Dirichlet and Drude conditions 
is more complex. Since monopoles I — occur only in 
the sector for m = 0, it is sufficient to reconsider the part 
of the energy coming from m = modes, and to trans- 
fer the Dirichlet energy unchanged for m =^ 0, giving the 
classical Casimir energy for the Drude case 



£(Dr) = E 



(Dr) 
m=0 



k R T 



i=i 



2Zln(l - Z 2l+1 ) , (5) 



where E^I_ I is given by Eq. 



2\ with N from Eq. \A\ with 



which we denote by No in the following. To compute the 
determinant in Eq. ([2]), we perform two transformations 
on the matrix that represents A*o : (i) A translation from 
the spherical wave basis Sc with origin C to a spherical 
wave basis Sf with origin F, see Fig.[l] (ii) A conversion 
from the basis Sf to the bispherical wave basis BSf with 
the same origin F. The translation (i) over the distance 
ZR corresponds to a similarity transformation with a 
matrix which has non-zero elements only for I > /' > 1, 
given by 

sg^s F{ {ZR) 1 - 1 ' I! [WTT 



(l-l')l 

The inverse of this matrix is obtained by ZR 
After this transformation, we can write A^ as 



N 



if / > i' 

Why. - (-I)') Z W ' +1 if I < V 



-ZR. 



(6) 



The conversion (ii) corresponds to a similarity transfor- 
mation with a matrix which has non-zero elements only 
for V > I > 1, given by 



(-1)' 

\/2T+T 



1+1/2 



n 



i\(i'-iy. 



With this transformation, the equivalence class can be 
expressed as 



No 



S w + (1 - Z 2 ){1 - Z 21 ) 



(7) 



m = and restricted to I, I' > 1, the equivalence class of 



The latter expression allows for a direct computation of 
det(l — No)- The first part in Eq. ^ is diagonal and 
can be easily factorized so that it yields a contribution 
(fe s T/2)E^iln(l - Z 2l+1 ) to the energy E^l- The 
second part in Eq. Q depends only the column index 
V and hence is a matrix with equal rows. For a matrix 
A of this type one has det(l — A) = 1 — trA Hence the 
contribution of the second part in Eq. ^ to the energy is 
given by the trace of this part divided by 1 — Z 2l+1 due to 
the factorization of the first part of 1 — No- Combining 
the two parts from Eq. Q with Eq. (J5j) we obtain the 
following exact expression for the classical Casimir energy 
of a Drude sphere and plate 



(Dr) ^bT 



£(Dr) = 



^(2Z + 1) ln(l - Z 2l+1 ) + In 1 - (1 - Z 2 ) £ Z 2l+1 - 



1=1 



, l -Z 21 

Z 2l +! 



(8) 



Applying similar arguments as in the Dirichlet case, it can again be shown that l?( Dr ) is an analytic function for 



4 



\Z\ < 1. We note that the effect of eliminating monopole 
fluctuations from the energy for Dirichlet conditions is 
not only a sum starting at I = 1 in E^ but the occur- 
rence of a second logarithmic term. This is a consequence 
of the coupling of monopolcs to higher order multipoles. 

Expressions at short distances — With exact expres- 
sions for the Casimir energies in Dirichlet and Drudc 
model available, one can compute explicitly the interac- 



tion in the experimentally important limit of short dis- 
tances L < R. Since this limit corresponds to Z close 
to unity, we compute the series in Eqs. ([3]), Q using the 
Abcl-Plana formula [2J. We set Z = cxp(— At) and expand 
for small n where /i = ln(l + £ + y/£(2 + £)), £ = L/R, in 
the sphere-plate geometry. For Dirichlet conditions, we 
obtain 



knT 



C(3)l 
2 



1 , 1 

— m u -\ — 
12 ^ 8 



7o 



31 



2880' 



725760 



(9) 



with the constant 70 = 0.174897 that is given by an integral [22], and for the Drude model we get 



e (Dt) = E (D) 



ln(7i — ln/i) + 



1 



-72 



ln/i 2 



1 73 - 74 In j" + In A 4 4 



6 —71 + In /i 



180 (-71+ In At) 2 



A 4 + <W 



(10) 



where the constants 71 = 1.270362, 72 = 1.35369, 
73 = 1.59409, 74 = 2.51153 are given by integrals that 
can be computed easily numerically. We used At as the 
variable for the expansion since it provides a very accu- 
rate result at even larger distances, see Fig. [2] A general 
feature of both the Dirichlet and Drude energy is their de- 
pendence on only In At and even powers of fx. This shows 
that the energies depend only on In £ and integer powers 
of £. When the force for the Dirichlet case is expanded 
in £, it becomes a Laurent series starting with l/£ 2 , i.e., 
there are no logarithmic terms for the force. For the 
Drude model, there are logarithmic terms in the force, 
and the most convenient form to express the short dis- 



tance expansion appears to be the one in Eq. (10). The 
leading correction to the PFA is the same term ~ In At in 
both models. However, at realistic distances, the double 



logarithmic term in Eq. ( 10 ) dominates and therefore the 



two models show rather distinct behavior, see Fig. [2] 

It is instructive to discuss the different behavior of 
the classical Casimir energies for Dirichlet and Drudc 
boundary conditions. For large separations, the interac- 
tions display different asymptotic behaviors for the two 
boundary conditions. In the experimentally most rele- 
vant sphere-plate geometry Eq. ^ predicts the known 
l/£ fall-off rate for the Dirichlet energy, while in the 
Drude case Eq. ^ gives the characteristic l/£ 3 decay 
|23j . The slower decay rate in the Dirichlet case re- 
sults from monopole contributions that are absent in 
the Drude case. While a Dirichlet scalar field is in gen- 
eral not expected to describe the Casimir interaction be- 
tween metals, in the high-temperature limit the differ- 
ence between the two universal interactions for Dirich- 
let and Drude conditions suggests an interesting phys- 
ical interpretation. In fact, the high-temperature limit 



of the Casimir interaction provides the dominant con- 
tribution to the full quantum Casimir interaction at fi- 
nite temperatures in the limit of sufficiently large sep- 
arations [HE]. To understand if the Dirichlet (l/£) 
or Drude (l/^ 3 ) decay describes experiments with con- 
ductors at large distances, one must bear in mind that 
the conductors used in Casimir force measurements are 
always connected to a charge reservoir to compensate 
stray charges that might otherwise be present on the sur- 
faces. In the static limit, the surface electric potential of 
such a conductor is constant (Dirichlet boundary condi- 
tions) and their total charge can fluctuate. In contrast, 
Drude metallic boundary conditions instead describe un- 
grounded charge-neutral conductors. From this we con- 
clude that the quantum Casimir force between grounded 
conductors at finite temperatures should decay according 
to the Dirichlet case at large distances since the energy 
is than dominated by the lowest Matsubara mode, i.e., 
the classical energy [HE]. To discriminate between the 
asymptotic l/£ and l/£ 3 decay, it is necessary to consider 
sphere-plate separations comparable to the sphere-radius 
or larger. Interestingly, our short distance expansions 
in the classical limit suggest a distinct behavior due to 
grounding at even shorter separations which, however, is 
certainly modified due to quantum corrections. 
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